Effective masses of diquarks 

P. Maris* 

Dept. of Physics, North Carolina State University, Raleigh, NC 27695-8202 

February 8, 2008 



Abstract 

We study meson and diquark bound states using the rainbow-ladder 
truncation of QCD's Dyson-Schwinger equations. The infrared strength of 
the rainbow-ladder kernel is described by two parameters. The ultraviolet 
behavior is fixed by the one-loop renormalization group behavior of QCD, 
which ensures the correct asymptotic behavior of the Bethe-Salpeter am- 
plitudes and brings important qualitative benefits. The diquark with the 
lowest mass is the scalar, followed by the axialvector and pseudoscalar 
diquark. This ordering can be anticipated from the meson sector. 



1 Introduction 

Mesons are color-singlet bound states of a quark and an antiquark. In addition 
to quark-antiquark bound states, one could also ask the question whether or 
not there are quark-quark bound states, by studying the corresponding Bethe- 
Salpeter equation [BSE] for bound states. Such states are of course not colorless 
in QCD, and are therefore expected to be confined, if they exist at all. Nev- 
ertheless, the masses of these states can serve as an indication for the relevant 
mass scales of quark-quark correlations. Such diquark correlations could play a 
role inside baryons: two quarks bound in a color antitriplet configuration can 
couple with a quark to form a color-singlet baryon. Indeed, recent studies of 
baryons as bound states of a quark and a (confined) diquark are quite successful 
in describing baryons [0, ||, ||, U . Also some lattice simulations || indicate the 
existence of correlations in diquark channels. 

Here we report ground state meson and color-antitriplet diquark masses 
based on the rainbow-ladder truncation of the set of Dyson-Schwinger equations 
[DSEs] || This covariant approach accommodates quark confinement and 
implements QCD's one-loop renormalization group behavior. The model we ap- 
ply has previously been shown to give an efficient description of the masses and 
electroweak decay constants of the light pseudoscalar and vector mesons fi ^| 
and their interactions |ll], 111], [ill [u| . This calculation will provide information 
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to guide the improvement of quark-diquark models of baryons and assists in 
developing an intuitive understanding of lattice simulations. 

2 Dyson— Schwinger Equations 

The DSE for the renormalized quark propagator in Euclidean space is 

r A 

Sip)' 1 = iZ 2 p' + Z 4 m q (ri + Z 1 g 2 D^(k)^S(q)ri(q : p) , (1) 

Jq 

where D^k = p — q) is the dressed-gluon propagator, F l v (q,p) the dressed- 
quark-gluon vertex with color-octet index i = 1, ...,8, and Zi and Z/± are 
the quark wave-function and mass renormalization constants. The notation 
= J d 4 q/ (2ir) 4 stands for a translationally invariant regularization of the 
integral, with A being the regularization mass-scale. The regularization is re- 
moved at the end of all calculations, by taking the limit A — > oo. We use the 
Euclidean metric where {7^, 7„} = 28^, 7^ = 7^ and a ■ b — J2i=i a ibi- The 
most general solution of Eq. ([!]) has the form S(p)- 1 = i^>A(p 2 )+ B(p 2 ) and is 
renormalized at spacelike p 2 according to A(^l 2 ) = 1 and B(p 2 ) = m q (n) with 
m q (p) being the current quark mass. 

2.1 Meson Bethe-Salpeter Equation 

Mesons are color-singlet quark-antiquark bound states. They are solutions of 
the homogeneous Bethe-Salpeter equation [BSE] for qq bound states 

TM(p+,P-)$ab = / K"^(p+,q+;q-,p-) 

Jq 

S^'(q + )Tlf(q + ,q_)6 cd S S ' 5 (q-), (2) 

where p+ = p + r]P and p- = p — (1 — 77) P are the outgoing and incoming 
quark momenta respectively, and q± is defined similarly, see also Fig. [l]. The 
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Figure 1: Meson Bethe-Salpeter Equation. 

greek superscripts are spinor indices, and the roman subscripts are color indices 
running from 1 to 3; for simplicity we consider one flavor here. The kernel 
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K is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. 

This equation has solutions at discrete values of P 2 = —m 2 M , where tum 
is the meson mass. Together with the canonical normalization condition for qq 
bound states, it completely determines Tm, the bound state Bethe-Salpeter am- 
plitude [BSA]. The different types of mesons, such as pseudo-scalar, vector, etc. 
are characterized by different Dirac structures. The most general decomposition 
for pseudoscalar bound states is 

r PS {q+,q-) = J 5 [lE(q 2 ,q-P;? 1 ) + P F(q 2 , q ■ P; V ) 

+ ft G{q 2 ,q ■ P; 77) + a MV P^q v H(q 2 , q ■ P; 77)] , (3) 

where the invariant amplitudes E, F, G and H are Lorentz scalar functions of q 2 
and q ■ P. Note also that these amplitudes explicitly depend on the momentum 
partitioning parameter r\. However, so long as Poincare invariance is respected, 
the resulting physical observables are independent of this parameter [§L [l^] . If 
77 = i , these amplitudes are appropriately odd or even in the charge parity odd 
quantity q ■ P for charge eigenstates. In the case of the h pion, for example, 
the amplitude G is odd in q ■ P, the others are even. 

The most general decomposition for scalar mesons can be obtained from 
that for a pseudoscalar mesons by dropping the 75. For the ++ , the amplitude 
F is odd in q ■ P, the others are even, if we use the same notation as for the 
pseudoscalars. A vector meson has more Dirac structures: a massive vector 
meson being transverse, the general decomposition of a such a BSA requires 
eight covariants ||, the dominant structure being 

K° m (q+,q-) = (s ii u-^-y^V 1 (q 2 ,q-P;r 1 ). (4) 

For the 1~~ p meson, the function Vi(q 2 , q ■ P; rj = i) is even in q ■ P. 

2.2 Diquark bound states 

Diquarks are quark-quark correlations. In QCD, with N c — 3, these states are 
necessarily colored and therefore believed to be confined. Two quarks can be 
coupled in either a color sextet or a color antitriplet. Single gluon exchange leads 
to an (effective) interaction that is attractive for diquarks in a color antitriplet 
configuration, but the interaction is repulsive in the color sextet channel JHJ [l5| . 
Furthermore, it is the diquark in a color antitriplet that can couple with a quark 
to form a color-singlet baryon. Thus we will only consider in diquarks in a color 
antitriplet configuration. 

Using the antisymmetric tensor e a & c as a representation of the antitriplet ]Tr| , 
the corresponding diquark bound states can be described by solutions of the 
homogeneous BSE 

,A 

^{p+,-P-Kbc = / K^2^(p+,<l+\-P-,-q-) 

Jq 

S^'(q+) TX?\q + -q-) e dec S 55 ' , (5) 
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where p+ = p + rjP and — p_ = — p + (1 — rj)P are now both outgoing quark 
momenta (and similarly for q + and — This is schematically depicted in 
Fig. |^; note the difference in the argument and the order of the indices of the 
second quark propagator and of the kernel K compared to the meson BSE, 
Eq. (0). With the help of the charge conjugation matrix C, which satisfies 
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Figure 2: Diquark Bethe-Salpeter equation. 
C 2 = —1 and C = (7 M ) T , we can rewrite the BSE for diquarks as juj 

(r D (p + , -p-)c) aP e abc = -e dec f K%} b e 5 (p + ,q +] -p., -<?_) C 5 " s C^' 



S^(q + )(r D (q + ,-q_)C)"' S %_) 



s5'6" 



(6) 



This equation resembles the meson BSE; in fact, using a skeleton expansion 
of the kernel K one can show that T£>(q +1 — q-)C satisfies a BSE whose Dirac 
structure is identical to that of the meson BSE for r m{<1+ ,9-)- Only its color 
structure is different, and that has important consequences |l7j ]. 

As in the case of mesons, the different types of diquarks are characterized 
by different Dirac structures. Since the intrinsic parity of a quark-quark pair 
is opposite to that of a quark-antiquark pair, a scalar diquark BSA, or rather 
r_D lS cai(9+, — Q—)C, has exactly the same form as a pseudoscalar meson BSA, 
Eq. (^). Similarly, an axialvector diquark BSA has the same decomposition as 
a vector meson BSA, and a pseudoscalar diquark BSA has the same form as a 
scalar meson BSA. 



2.3 Two-body bound states in rainbow-ladder truncation 

For practical calculations, we utilize the rainbow truncation of the quark DSE 

Z 1 g 2 D flu (k)ri(q,p) - g(k 2 )D^ c (k) lu ^- 7 (7) 
in combination with the ladder truncation of the BSE 
\y+>H+i'i-ij j -) ~^ 

A 



T 



-Q(k z 



(7m) 7 D^{k) 



(7, 



(8) 
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where D ™ e (k = p — q = p± — q±) is the free gluon propagator in Landau gauge. 
These two truncations are consistent in the sense that the combination produces 
vector and axialvector vertices satisfying the respective Ward-Takahashi iden- 
tities Q . In the axial case, this ensures that in the chiral limit the ground state 
pseudoscalar mesons are the massless Goldstone bosons associated with chiral 
symmetry breaking; with nonzero current quark masses it leads to a generaliza- 
tion of the Gell-Man-Oaks-Renner relation C^l . In the vector case, this ensures 
conservation of the electromagnetic current |10| . 

If we insert this ladder kernel, Eq. (||), in the meson BSE, we get the ladder 
BSE for the meson BSA 

fA 



y y 



5(fc 2 ) D^{k) 

a/3 



/ N Up 

\l^S{q+)T M {q + ,q^) ) 7^ , (9) 
For the diquark bound states we obtain 

(T D (p + ,-p_)c) a \ abc = £dec ^jj d (j) b fa*?) 

( 7p S(q+) T D (q + , - g _) C %_) 7/1 ) ^ . (10) 

Comparing Eq. (|^) with Eq. ([l0|), we see that in ladder truncation the only 
difference between the meson BSE for Tm and BSE for a color-antitriplet di- 
quark, or rather, for TjjC, lies in the color factors. For the mesons, we have 
Tr c [4r 4j-] = 4, and thus we get a factor of —4/3 in front of the BSE integral. 
For the diquarks in a color antitriplet configuration on the other hand, we have 

y) od (y) be ed ~ - (ID 

Thus, in ladder truncation the effective interaction for the diquarks is reduced 
by a factor of two compared to the interaction in the meson channel. Beyond 
ladder truncation a more complex algebraic relation persists between the meson 
and diquark BSEs [fl7| . We therefore expect that the diquarks with the lowest 
mass are the scalar diquarks (the diquark partners of the pseudoscalar mesons), 
followed by the axialvector and pseudoscalar diquarks |jl4j| . 

The rainbow-ladder truncation is particularly suitable for the flavor octet 
pseudoscalar mesons and for the vector mesons, since the next-to-leading-order 
contributions in a quark-gluon skeleton graph expansion have a significant 
amount of cancellation between repulsive and attractive corrections [jl7| |l9| ]. 
However, the BSE for diquarks, and also the scalar meson BSE, receive large 
repulsive corrections from these next-to-leading-order contributions. This will 
significantly change the scalar meson mass and corresponding BSA; the diquark 
bound states found in ladder truncation will in fact disappear altogether from 
the spectrum |l7], [l9[ ^0|. Nevertheless, correlations will persist in the diquark 
channels and diquark masses found in ladder truncation can serve as a guide for 
parametrizing those correlations. 
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3 Numerical calculations 



3.1 Model interaction 



For our numerical calculations we employ a model for the effective coupling 
Q{k 2 ) that has been developed through an efficient description of the masses and 
decay constants of the light pseudoscalar and vector mesons This effective 
qq interaction is constrained by perturbative QCD in the ultraviolet and has a 
phenomenological infrared behavior. The form is 



Q(k 2 
k 2 



iw 2 Dk 2 



^ 2 lm T(k 2 



I In 



1 



k 2 /A 2 



QCD) 



(12) 



with 7 m = 12/(33 - 2N f ) and T(s) = (1 - exp j^)/s. The ultraviolet behav- 
ior is chosen to be that of the QCD running coupling a(k 2 ); the ladder-rainbow 
truncation then generates the correct perturbative QCD structure of the DSE- 
BSE system of equations. The first term implements the strong infrared support 
in the region < k 2 < 1 GeV 2 phenomenologically required |y] to produce a 
realistic value for the chiral condensate. We use m t — 0.5 GeV, r = e 2 — 1, 
Nf = 4, Aqcd = 0.234 GeV, and a renormalization scale /i = 19 GeV which 
is well into the perturbative domain ||, The remaining parameters, 

uj = 0.4 GeV and D — 0.93 GeV 2 along with the quark masses, are fitted to 
give a good description of the chiral condensate, ra^j K and f T . 

Within this model, the quark propagator reduces to the one-loop pertur- 
bative QCD propagator in the ultraviolet region. In particular, the dynamical 
mass function M(p 2 ) = B(p 2 )/A(p 2 ) behaves like 



M(p 2 



rh q 


o 


r 


r 


. QCD _ 





(13) 



where rh q is the renormalization-point-independent explicit chiral-symmetry- 
breaking mass. In the chiral limit the behavior is qualitatively different 



M cMra i(p 2 ) 



2tt 2 7i) 



(qq)° 



p 2 ( i In 



l-7n 



(qq) lt = -3Z 4 [ Tr[5 chiral (p)] 

J a 



is at one-loop level 



(14) 



with {qq}° the renormalization-point-independent chiral condensate Its re- 
lation to the renormalization-dependent condensate 



(15) 



(\nn/A QCD y™(qq)°. 



(16) 



G 



As demonstrated clearly in Ref. Eq], this behavior of the mass function in the 
chiral limit is a keystone of the microscopic realization of Goldstone's theorem 
in QCD: the quark condensate (and hence the pion mass) and other observables 
are materially dependent on the ultraviolet properties of the interaction. 

In the infrared region both the wave function renormalization Z(p 2 ) — 
1/A(p 2 ) and the dynamical mass function deviate significantly from the pertur- 
bative behavior, due to chiral symmetry breaking. Recent comparisons [p2| , p3| 
of results from this rainbow DSE model to lattice QCD simulations [£4|] pro- 
vide semiquantitative confirmation of the behavior generated by the present 
DSE model: a significant enhancement of M(p 2 ) and a material enhancement 
of A(p 2 ) below 1 GeV 2 . 

The vector meson masses and electroweak decay constants produced by this 
model are in good agreement with experiments The model's prediction 
for the pion charge form factor F n (Q 2 ) is confirmed by the most recent Jlab 
data [p5| . The kaon charge radii and electromagnetic form factors are also 
described well |L0|] , as is the weak Ki 3 decay The strong decays of the 

vector mesons into a pair of pseudoscalar mesons also agree reasonably well 
with experiments JlJ, [2(|. The performance of this model for deep inelastic 
scattering phenomena can be gauged from that of a simplified version that has 
recently produced reasonable results for the pion valence quark distribution 
amplitude [ p7[ . 



3.2 Ground state mass spectrum of qq and qq states 

With this model we have calculated the masses of the lightest two-body bound 
states in the u, d, s quark sector: pseudoscalar, scalar, and vector mesons, and 
scalar, pseudoscalar, and axial vector diquarks. As far as the flavor-singlet 
mesons are concerned, with the rainbow-ladder kernel only flavor-singlet states 
with "ideal mixing" are possible. This feature indicates a deficiency of the 
rainbow-ladder truncation for studying flavor-singlet pseudoscalar and scalar 
mesons (the flavor-singlet vector mesons are very close to ideal mixing). 

Regarding the flavor structure of the diquarks, the Pauli principle prohibits 
the existence of flavor-singlet color-antitriplet scalar and pseudoscalar diquarks; 
the only possible ss diquarks are axialvector diquarks j^, |l6). Nevertheless, 
in model studies such as this, one can consider the up/down quark masses to 
be equal to the strange quark mass, and calculate the corresponding meson 
and diquark masses. This indicates the dependence of the diquark masses on 
the current quark masses and allows for a comparison with other calculations. 
However, one should keep in mind that there are no scalar nor pseudoscalar ss 
diquarks in 3-flavor QCD with realistic quark masses. 

The bound state BSE has solutions at discrete values of P 2 = —to 2 , where 
m is the mass of the bound state. In order to determine the bound state masses, 
we introduce a parameter A(P 2 ), and turn the BSE into an eigenvalue problem 
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for A(P 2 ) 



A g(k 2 )D^ c (k)^S(q + )T(q + ,q^S(q^^ = A(P 2 ) T(p + ,p^) . (17) 

This equation has solutions for all P 2 , but in general these solutions do not 
correspond to bound state solutions of the original BSE, Eqs. @ and ( |l0|) . 
Meson bound states with mass to are described by solutions of Eq. ( p.7| ) 
r(p + ,p_) = Tm{p+,P-) with A(P 2 = — to 2 ) = 1, whereas diquark bound states 
are described by solutions T(p + ,p_) = Te>(p+,— P-)C with A(— to 2 ) = 2 (re- 
member that the color structure in the diquark BSE leads to a factor | instead 
of | , see Eq. In Fig. || we display the eigenvalues as function of P 2 in 

the up/down sector. From this figure we can simply read of the values of P 2 
that give us the mass of a bound state. 



• pseudoscalar meson / scalar diquark 
■ vector meson / axial-vector diquark 
+ scalar meson / pseudoscalar diquark 




-1.2 -0.8 -0.4 0.4 0.8 1.2 1.6 2 

P 2 [GeV 2 ] 



Figure 3: The eigenvalue A(P 2 ) as function of P 2 in the up/down sector: A = 1 
signals a meson bound state, whereas A = 2 corresponds to a diquark bound 
state. 

Since we started from a Euclidean formulation, we have to make an ana- 
lytic continuation to negative values of P 2 in order to determine the bound 
state solutions. For small masses, this is no problem; however, for larger masses 
this becomes numerically more and more cumbersome. Moreover, the quark 
propagator in rainbow truncation often has singularities at complex values of 
p 2 p8| , which could be interpreted as a signal for confinement ]29| ], Once the 
integration domain sampled in the BSE includes these singularities, one has to 
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specify how to deal with these singularities. Alternatively, one can extrapolate 
the eigenvalues A(P 2 ) from the region where the analytic continuation is unam- 
biguous to A = 1 for mesons, and to A = 2 for diquarks, whenever necessary. We 
followed the latter procedure, using several different extrapolating functions to 
reduce the numerical uncertainty. For the bound states with the largest mass, 
we estimate the error due to this extrapolation to be a few percent. 
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Table 1: Overview of our results for the meson and diquark masses for different current 
quark masses, all in GeV, compared with the results from a separable model terj and with 
quenched lattice QCD Note that the Pauli principle prohibits the existence of flavor- 
singlet color-antitriplet scalar and pseudoscalar diquarks; the results for the m q = m q r = m a 
scalar and pseudoscalar diquark masses are only included to indicate the dependece on the 
current quark mass and for comparison with lattice data. Our numerical errors are of the 
order of 1%, with the exception of the heavier-mass states, which have an extrapolation error 
as indicated. 

Our results are displayed in Table [I], and compared with results from 
rainbow-ladder truncation of the DSEs using a separable Ansatz for the inter- 
action [|o| and with lattice data where available. Note that our calculation 
is done at realistic quark masses, whereas the lattice results are extrapolated 
to the chiral limit. For the strange quark we have taken the lattice data with 
my = Tricf,. Compared to the lattice data, we find a larger splitting between 
the scalar and the axialvector diquark; also in the ss sector we find a signifi- 
cantly larger splitting between the pseudoscalar and vector meson. The lattice 
simulations do not report on scalar meson or pseudoscalar diquark masses. 

Our results are qualitatively the same as the results of the separable model 
of Ref. |3(J, with the exception of the scalar mesons and pseudoscalar diquarks. 
We find the scalar and vector meson masses to be within 10% of each other, not 
only in the u/d sector, but also in the us and ss sectors. Similarly, we also find 
the pseudoscalar diquark and axialvector diquark masses to be within about 
10% of each other, whereas in the separable model the pseudoscalar diquark is 
much heavier than the axialvector diquark, or not bound at all. Also notice 
that the scalar meson and pseudoscalar diquark masses increase more rapidly 
with the quark mass than the vector meson and scalar diquark masses: e.g. the 
m CT < m p but TOp+ > 

In Table | we also give the masses obtained using only the dominant co- 
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variant for the different bound states; e.g., retaining only E in Eq. (||). This 
approximation changes the mass by about 20% both for the mesons and for the 
diquarks. Using the dominant covariants only tends to increase the difference 
between pseudoscalar meson masses and vector and scalar masses; we see a 
similar increase in the difference in the diquark channels. 

The meson masses seem to be rather independent of the details of the ef- 
fective interaction, as long as the interaction generates the observed amount of 
chiral symmetry breaking, as can be seen from Table |. The parameters to and 
D were fitted in Ref. U to reproduce /„• and the chiral condensate. Within 
this parameter range, the vector meson masses are almost independent of the 
parameters, and the scalar meson masses change by only 10%. The diquark 
masses are more sensitive to details of the effective interaction. Nevertheless, 
the lightest diquark is the scalar diquark, independent of the parameters, with 
the axialvector and pseudoscalar diquarks being about 150 to 250 MeV heavier. 
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Table 2: Parameter dependence of the masses (in GeV) in the up/down sector. 



3.3 Bethe— Salpeter amplitudes 

The corresponding BSAs are similar for the different bound states. In Fig. [| we 
have displayed the leading Chebyshev moments 

f(q 2 ) = f sin 2 d0 f(q 2 , qP COS0; V = ±) (18) 
Jo 

of the amplitudes associated with the canonical Dirac structure for the various 
mesons and diquarks in the up/down sector, all normalized to /(0) = 1: E 
for pseudoscalar and scalar bound states, see Eq. (||), and Vi for vector and 
axialvector bound states, see Eq. (|^) . 

In the infrared region, < q 2 < 2 GeV 2 , the canonical covariants can 
reasonably well be approximated by exponentials f(q 2 ) — exp(—q 2 /uj 2 ) with 
0.74 GeV < lu < 1.01 GeV, depending on the meson or diquark, see Tablet! At 
large q 2 , they all fall off like c/q 2 up to calculable logarithmic corrections ||, 
It is thus apparent that in all three channels, the diquark BSAs are narrower 
in momentum space than the corresponding meson BSAs, and thus wider in 
coordinate space. This is what one would expect, since the diquarks are less 
bound. It is intriguing that the scalar meson BSA is broader in momentum 
space than the pion BSA, and thus narrower in coordinate space. One expects 
the physical a state to include significant pion cloud contributions and to thus be 
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Figure 4: Leading Chebyshcv moments of canonical amplitudes in the up/down 
sector. 



mesons diquarks 
state n p a 0+1+0" 

u [GcV] 0.84 0.79 1.01 0.76 0.74 0.82 



Table 3: Comparison of the BSAs for the up/down states, fitted to exp(— q 2 /u> 2 ) in the 
infrared region. 

significantly broader in coordinate space than the pion. In ladder BSE approx- 
imation, such physical contributions are not accommodated; also other terms 
beyond rainbow-ladder truncation are known to be important in the scalar chan- 
nel. Clearly, additional study beyond ladder truncation is required in the scalar 
channel. 

4 Concluding remarks 

We have calculated the meson and color-antitriplet diquark masses using DSEs 
in rainbow-ladder truncation. The model parameters have previously been fixed 
by the calculation of m ff , rajf, and the chiral condensate, and it leads 
to a good description of the light pscudoscalar and vector mesons and their 
interactions. Using this model, we find that the lightest diquark bound states 
in a color antitriplet configuration is a scalar diquark, with a mass of about 



11 



800 MeV in the case of up and down quarks. The axialvector and pseudoscalar 
diquarks are about 150 to 250 MeV heavier than the scalar diquark. In the 
meson channel, we find a scalar meson with a mass slightly below the vector 
meson for up/down quarks, but slightly above the vector meson mass for strange 
quarks. These results are consistent with the notion of a (broad) cr-meson with 
a mass around 600 MeV, and with the mass splitting between the scalar and 
vector mesons in heavy quarkonia. The corresponding BSAs are all similar in 
shape, and can in the infrared region be described by a gaussian. For large 
relative momenta, they fall off like 1/p 2 , with calculable logarithmic corrections 
which is characteristic of QCD. 

The ladder truncation is particularly suitable for the pseudoscalar mesons 
and for the vector mesons, because of cancellations in the higher-order contribu- 
tions to the kernel. However, diquarks, and also the scalar meson, are expected 
to receive significant corrections from higher order terms in the kernel. This 
can significantly change the bound state mass and corresponding BSA; in par- 
ticular the diquarks will no longer be a true bound state. Thus one should 
only consider these masses as an indication for the relevant mass scales of the 
diquark correlations. As such they should be useful in e.g. studies of baryons 
as a quark-diquark system. 

Typically, studies of baryons as a quark-diquark bound state include both 
scalar and axialvector diquarks in a color antitriplet, but no pseudoscalar di- 
quarks. The present study indicates that the mass scale of pseudoscalar diquarks 
is similar to that of axialvector diquarks. It may therefore be interesting to ex- 
plore the influence of pseudoscalar correlations on quark-diquark models for 
baryons such as those of Refs. [Q, |[ |[ ||. 
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